New insights into the molecular basis of spinal neurofibromatosis type 1.

Spinal neurofibromatosis (SNF) is a form of neurofibromatosis type 1 (NF1) characterized by bilateral neurofibromas involving all spinal roots. The pathogenic mechanisms determining the SNF form are currently unknown. To verify the presence of genetic variants possibly related to SNF or classic NF1, we studied 106 sporadic NF1 and 75 SNF patients using an NGS panel of 286 genes encoding RAS pathway effectors and neurofibromin interactors and evaluated the expression of syndecans (SDC1, SDC2, SDC3, SDC4), the NF1 3' tertile interactors, by quantitative real-time PCR. We previously identified 75 and 106 NF1 variants in SNF and NF1 cohorts, respectively. The analysis of the distribution of pathogenic NF1 variants in the three NF1 tertiles showed a significantly higher prevalence of NF1 3' tertile mutations in SNF than in the NF1 cohort. We hypothesized a potential pathogenic significance of the 3' tertile NF1 variants in SNF. The analysis of syndecan expression on PBMCs RNAs from 16 SNF, 16 classic NF1 patients and 16 healthy controls showed that the expression levels of SDC2 and SDC3 were higher in SNF and NF1 patients than in controls; moreover, SDC2, SDC3 and SDC4 were significantly over expressed in patients mutated in the 3' tertile compared to controls. Two different mutational NF1 spectra seem to characterize SNF and classic NF1, suggesting a pathogenic role of NF1 3' tertile and its interactors, syndecans, in SNF. Our study, providing new insights on a possible role of neurofibromin C-terminal in SNF, could address effective personalized patient management and treatments.

A multi-step genomic approach prioritized TBKBP1 gene as relevant for multiple sclerosis susceptibility.

BACKGROUND: Over 200 genetic loci have been associated with multiple sclerosis (MS) explaining ~ 50% of its heritability, suggesting that additional mechanisms may account for the "missing heritability" phenomenon. OBJECTIVE: To analyze a large cohort of Italian individuals to identify markers associated with MS with potential functional impact in the disease. METHODS: We studied 2571 MS and 3234 healthy controls (HC) of continental Italian origin. Discovery phase included a genome wide association study (1727 MS, 2258 HC), with SNPs selected according to their association in the Italian cohort only or in a meta-analysis of signals with a cohort of European ancestry (4088 MS, 7144 HC). Top associated loci were then tested in two Italian cohorts through array-based genotyping (903 MS, 884 HC) and pool-based target sequencing (588 MS, 408 HC). Finally, functional prioritization through conditional eQTL and mQTL has been performed. RESULTS: Top associated signals overlap with already known MS loci on chromosomes 3 and 17. Three SNPs (rs4267364, rs8070463, rs67919208), all involved in the regulation of TBKBP1, were prioritized to be functionally relevant. CONCLUSIONS: No evidence of novel signal of association with MS specific for the Italian continental population has been found; nevertheless, two MS loci seems to play a relevant role, raising the interest to further investigations for TBKBP1 gene.

A Translational Approach to Spinal Neurofibromatosis: Clinical and Molecular Insights from a Wide Italian Cohort.

Spinal neurofibromatosis (SNF), a phenotypic subclass of neurofibromatosis 1 (NF1), is characterized by bilateral neurofibromas involving all spinal roots. In order to deepen the understanding of SNF's clinical and genetic features, we identified 81 patients with SNF, 55 from unrelated families, and 26 belonging to 19 families with at least 1 member affected by SNF, and 106 NF1 patients aged >30 years without spinal tumors. A comprehensive NF1 mutation screening was performed using NGS panels, including NF1 and several RAS pathway genes. The main features of the SNF subjects were a higher number of internal neurofibromas (p < 0.001), nerve root swelling (p < 0.001), and subcutaneous neurofibromas (p = 0.03), while hyperpigmentation signs were significantly less frequent compared with the classical NF1-affected cohorts (p = 0.012). Fifteen patients underwent neurosurgical intervention. The histological findings revealed neurofibromas in 13 patients and ganglioneuromas in 2 patients. Phenotypic variability within SNF families was observed. The proportion of missense mutations was higher in the SNF cases than in the classical NF1 group (21.40% vs. 7.5%, p = 0.007), conferring an odds ratio (OR) of 3.34 (CI = 1.33−10.78). Two unrelated familial SNF cases harbored in trans double NF1 mutations that seemed to have a subclinical worsening effect on the clinical phenotype. Our study, with the largest series of SNF patients reported to date, better defines the clinical and genetic features of SNF, which could improve the management and genetic counseling of NF1.

Genomic and functional evaluation of TNFSF14 in multiple sclerosis susceptibility.

Among multiple sclerosis (MS) susceptibility genes, the strongest non-human leukocyte antigen (HLA) signal in the Italian population maps to the TNFSF14 gene encoding LIGHT, a glycoprotein involved in dendritic cell (DC) maturation. Through fine-mapping in a large Italian dataset (4,198 patients with MS and 3,903 controls), we show that the TNFSF14 intronic SNP rs1077667 is the primarily MS-associated variant in the region. Expression quantitative trait locus (eQTL) analysis indicates that the MS risk allele is significantly associated with reduced TNFSF14 messenger RNA levels in blood cells, which is consistent with the allelic imbalance in RNA-Seq reads (P < 0.0001). The MS risk allele is associated with reduced levels of TNFSF14 gene expression (P < 0.01) in blood cells from 84 Italian patients with MS and 80 healthy controls (HCs). Interestingly, patients with MS are lower expressors of TNFSF14 compared to HC (P < 0.007). Individuals homozygous for the MS risk allele display an increased percentage of LIGHT-positive peripheral blood myeloid DCs (CD11c+, P = 0.035) in 37 HCs, as well as in in vitro monocyte-derived DCs from 22 HCs (P = 0.04). Our findings suggest that the intronic variant rs1077667 alters the expression of TNFSF14 in immune cells, which may play a role in MS pathogenesis.

Contribution of Rare and Low-Frequency Variants to Multiple Sclerosis Susceptibility in the Italian Continental Population.

Genome-wide association studies identified over 200 risk loci for multiple sclerosis (MS) focusing on common variants, which account for about 50% of disease heritability. The goal of this study was to investigate whether low-frequency and rare functional variants, located in MS-established associated loci, may contribute to disease risk in a relatively homogeneous population, testing their cumulative effect (burden) with gene-wise tests. We sequenced 98 genes in 588 Italian patients with MS and 408 matched healthy controls (HCs). Variants were selected using different filtering criteria based on allelic frequency and in silico functional impacts. Genes showing a significant burden (n = 17) were sequenced in an independent cohort of 504 MS and 504 HC. The highest signal in both cohorts was observed for the disruptive variants (stop-gain, stop-loss, or splicing variants) located in EFCAB13, a gene coding for a protein of an unknown function (p < 10-4). Among these variants, the minor allele of a stop-gain variant showed a significantly higher frequency in MS versus HC in both sequenced cohorts (p = 0.0093 and p = 0.025), confirmed by a meta-analysis on a third independent cohort of 1298 MS and 1430 HC (p = 0.001) assayed with an SNP array. Real-time PCR on 14 heterozygous individuals for this variant did not evidence the presence of the stop-gain allele, suggesting a transcript degradation by non-sense mediated decay, supported by the evidence that the carriers of the stop-gain variant had a lower expression of this gene (p = 0.0184). In conclusion, we identified a novel low-frequency functional variant associated with MS susceptibility, suggesting the possible role of rare/low-frequency variants in MS as reported for other complex diseases.

A reliable strategy for single-cell RNA sequencing analysis using cryoconserved primary cortical cells.

BACKGROUND: The application of single-cell RNA sequencing (scRNASeq) represents a unique approach to identify hundreds to millions of cells in mammalian cortical multilayers at different stages of embryogenesis. ScRNASeq technology applied to neurological studies requires the use of fresh starting materials because standard cryopreservation methods do not guarantee high viability of cortical primary cells derived from dissected brain areas. NEW METHOD: Here we set up and validate an innovative strategy to perform scRNASeq studies in cryopreserved primary cortical cells isolated from E15.5 mouse embryo. In order to freeze cortical primary cells, we have employed Neurostore, a medium able to guarantee high viability and cell composition of embryonic cortex after thawing. COMPARISON WITH EXISTING METHODS: We showed for the first time the possibility to run scRNASeq experiments on primary cortical cells in an off-line set-up, ensuring cellular integrity and diversity. RESULTS: By trypan blue assay and flow cytometry analysis, we found that Neurostore-cryopreserved cortical cells showed approximately 95 % of viability. Satisfactory RNA recovery and cDNA libraries were achieved. Transcriptome sequencing of 35,763 cryoconserved single cells yielded a robust data-set, identifying 25 cell clusters in three biological samples. Prevalence of peculiar neural populations before and after the cryopreservation-resuscitation procedure was verified by marker gene expression and immunofluorescence analysis. CONCLUSIONS: Our findings support the evidence that frozen primary cortical cells can be successfully employed in scRNASeq experiments allowing an unprecedented flexibility in experimental procedures, such as sample preparation and subsequent processing steps performed in different locations.

The Genome-Wide Impact of Nipblb Loss-of-Function on Zebrafish Gene Expression.

Transcriptional changes normally occur during development but also underlie differences between healthy and pathological conditions. Transcription factors or chromatin modifiers are involved in orchestrating gene activity, such as the cohesin genes and their regulator NIPBL. In our previous studies, using a zebrafish model for nipblb knockdown, we described the effect of nipblb loss-of-function in specific contexts, such as central nervous system development and hematopoiesis. However, the genome-wide transcriptional impact of nipblb loss-of-function in zebrafish embryos at diverse developmental stages remains under investigation. By RNA-seq analyses in zebrafish embryos at 24 h post-fertilization, we examined genome-wide effects of nipblb knockdown on transcriptional programs. Differential gene expression analysis revealed that nipblb loss-of-function has an impact on gene expression at 24 h post fertilization, mainly resulting in gene inactivation. A similar transcriptional effect has also been reported in other organisms, supporting the use of zebrafish as a model to understand the role of Nipbl in gene regulation during early vertebrate development. Moreover, we unraveled a connection between nipblb-dependent differential expression and gene expression patterns of hematological cell populations and AML subtypes, enforcing our previous evidence on the involvement of NIPBL-related transcriptional dysregulation in hematological malignancies.

Digenic inheritance of subclinical variants in Noonan Syndrome patients: an alternative pathogenic model?

Noonan syndrome (NS) is an autosomal-dominant disorder with variable expressivity and locus heterogeneity. Despite several RAS pathway genes were implicated in NS, 20-30% of patients remain without molecular diagnosis, suggesting the involvement of further genes or multiple mechanisms. Eight patients out of 60, negative for conventional NS mutation analysis, with heterogeneous NS phenotype were investigated by means of target resequencing of 26 RAS/MAPK pathway genes. A trio was further characterized by means of whole-exome sequencing. Protein modeling and in silico prediction of protein stability allowed to identify possible pathogenic RAS pathway variants in four NS patients. A new c.355T>C variant in LZTR1 was found in patient 43. Two patients co-inherited variants in LRP1 and LZTR1 (patient 53), or LRP1 and SOS1 genes (patient 67). The forth patient (56) carried a compound heterozygote of RASAL3 gene variants and also an A2ML1 variant. While these subclinical variants are singularly present in healthy parents, they co-segregate in patients, suggesting their addictive effect and supporting a digenic inheritance, as alternative model to a more common monogenic transmission. The ERK1/2 and SAPK/JNK activation state, assessed on immortalized lymphocytes from patients 53 and 67 showed highest phosphorylation levels compared to their asymptomatic parents. These findings together with the lack of their co-occurrence in the 1000Genomes database strengthen the hypothesis of digenic inheritance in a subset of NS patients. This study suggests caution in the exclusion of subclinical variants that might play a pathogenic role providing new insights for alternative hereditary mechanisms.

Author Correction: Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A novel homozygous mutation in the TRDN gene causes a severe form of pediatric malignant ventricular arrhythmia.

BACKGROUND: Triadin is a protein expressed in cardiac and skeletal muscle that has an essential role in the structure and functional regulation of calcium release units and excitation-contraction coupling. Mutations in the triadin gene (TRDN) have been described in different forms of human arrhythmia syndromes with early onset and severe arrhythmogenic phenotype, including triadin knockout syndrome. OBJECTIVE: The purpose of this study was to characterize the pathogenetic mechanism underlying a case of severe pediatric malignant arrhythmia associated with a defect in the TRDN gene. METHODS: We used a trio whole exome sequencing approach to identify the genetic defect in a 2-year-old boy who had been resuscitated from sudden cardiac arrest and had frequent episodes of ventricular fibrillation and a family history positive for sudden death. We then performed in vitro functional analysis to investigate possible pathogenic mechanisms underlying this severe phenotype. RESULTS: We identified a novel homozygous missense variant (p.L56P) in the TRDN gene in the proband that was inherited from the heterozygous unaffected parents. Expression of a green fluorescent protein (GFP)-tagged mutant human cardiac triadin isoform (TRISK32-L56P-GFP) in heterologous systems revealed that the mutation alters protein dynamics. Furthermore, when co-expressed with the type 2 ryanodine receptor, caffeine-induced calcium release from TRISK32-L56P-GFP was relatively lower compared to that observed with the wild-type construct. CONCLUSION: The results of this study allowed us to hypothesize a pathogenic mechanism underlying this rare arrhythmogenic recessive form, suggesting that the mutant protein potentially can trigger arrhythmias by altering calcium homeostasis.

Genome sequencing of Prototheca zopfii genotypes 1 and 2 provides evidence of a severe reduction in organellar genomes.

Prototheca zopfii (P. zopfii, class Trebouxiophyceae, order Chlorellales, family Chlorellaceae), a non-photosynthetic predominantly free-living unicellular alga, is one of the few pathogens belonging to the plant kingdom. This alga can affect many vertebrate hosts, sustaining systemic infections and diseases such as mastitis in cows. The aim of our work was to sequence and assemble the P. zopfii genotype 1 and genotype 2 mitochondrial and plastid genomes. Remarkably, the P. zopfii mitochondrial (38 Kb) and plastid (28 Kb) genomes are models of compaction and the smallest known in the Trebouxiophyceae. As expected, the P. zopfii genotype 1 and 2 plastid genomes lack all the genes involved in photosynthesis, but, surprisingly, they also lack those coding for RNA polymerases. Our results showed that plastid genes are actively transcribed in P. zopfii, which suggests that the missing RNA polymerases are substituted by nuclear-encoded paralogs. The simplified architecture and highly-reduced gene complement of the P. zopfii mitochondrial and plastid genomes are closer to those of P. stagnora and the achlorophyllous obligate parasite Helicosporidium than to those of P. wickerhamii or P. cutis. This similarity is also supported by maximum likelihood phylogenetic analyses inferences. Overall, the P. zopfii sequences reported here, which include nuclear genome drafts for both genotypes, will help provide both a deeper understanding of the evolution of Prototheca spp. and insights into the corresponding host/pathogen interactions.

A Census of Tandemly Repeated Polymorphic Loci in Genic Regions Through the Comparative Integration of Human Genome Assemblies.

Polymorphic Tandem Repeat (PTR) is a common form of polymorphism in the human genome. A PTR consists in a variation found in an individual (or in a population) of the number of repeating units of a Tandem Repeat (TR) locus of the genome with respect to the reference genome. Several phenotypic traits and diseases have been discovered to be strongly associated with or caused by specific PTR loci. PTR are further distinguished in two main classes: Short Tandem Repeats (STR) when the repeating unit has size up to 6 base pairs, and Variable Number Tandem Repeats (VNTR) for repeating units of size above 6 base pairs. As larger and larger populations are screened via high throughput sequencing projects, it becomes technically feasible and desirable to explore the association between PTR and a panoply of such traits and conditions. In order to facilitate these studies, we have devised a method for compiling catalogs of PTR from assembled genomes, and we have produced a catalog of PTR for genic regions (exons, introns, UTR and adjacent regions) of the human genome (GRCh38). We applied four different TR discovery software tools to uncover in the first phase 55,223,485 TR (after duplicate removal) in GRCh38, of which 373,173 were determined to be PTR in the second phase by comparison with five assembled human genomes. Of these, 263,266 are not included by state-of-the-art PTR catalogs. The new methodology is mainly based on a hierarchical and systematic application of alignment-based sequence comparisons to identify and measure the polymorphism of TR. While previous catalogs focus on the class of STR of small total size, we remove any size restrictions, aiming at the more general class of PTR, and we also target fuzzy TR by using specific detection tools. Similarly to other previous catalogs of human polymorphic loci, we focus our catalog toward applications in the discovery of disease-associated loci. Validation by cross-referencing with existing catalogs on common clinically-relevant loci shows good concordance. Overall, this proposed census of human PTR in genic regions is a shared resource (web accessible), complementary to existing catalogs, facilitating future genome-wide studies involving PTR.

T cell neoepitope discovery in colorectal cancer by high throughput profiling of somatic mutations in expressed genes.

OBJECTIVE: Patient-specific (unique) tumour antigens, encoded by somatically mutated cancer genes, generate neoepitopes that are implicated in the induction of tumour-controlling T cell responses. Recent advancements in massive DNA sequencing combined with robust T cell epitope predictions have allowed their systematic identification in several malignancies. DESIGN: We undertook the identification of unique neoepitopes in colorectal cancers (CRCs) by using high-throughput sequencing of cDNAs expressed by standard cancer cell cultures, and by related cancer stem/initiating cells (CSCs) cultures, coupled with a reverse immunology approach not requiring human leukocyte antigen (HLA) allele-specific epitope predictions. RESULTS: Several unique mutated antigens of CRC, shared by standard cancer and related CSC cultures, were identified by this strategy. CD8+ and CD4+ T cells, either autologous to the patient or derived from HLA-matched healthy donors, were readily expanded in vitro by peptides spanning different cancer mutations and specifically recognised differentiated cancer cells and CSC cultures, expressing the mutations. Neoepitope-specific CD8+ T cell frequency was also increased in a patient, compared with healthy donors, supporting the occurrence of clonal expansion in vivo. CONCLUSIONS: These results provide a proof-of-concept approach for the identification of unique neoepitopes that are immunogenic in patients with CRC and can also target T cells against the most aggressive CSC component.

Next Generation Sequencing of Pooled Samples: Guideline for Variants' Filtering.

Sequencing large number of individuals, which is often needed for population genetics studies, is still economically challenging despite falling costs of Next Generation Sequencing (NGS). Pool-seq is an alternative cost- and time-effective option in which DNA from several individuals is pooled for sequencing. However, pooling of DNA creates new problems and challenges for accurate variant call and allele frequency (AF) estimation. In particular, sequencing errors confound with the alleles present at low frequency in the pools possibly giving rise to false positive variants. We sequenced 996 individuals in 83 pools (12 individuals/pool) in a targeted re-sequencing experiment. We show that Pool-seq AFs are robust and reliable by comparing them with public variant databases and in-house SNP-genotyping data of individual subjects of pools. Furthermore, we propose a simple filtering guideline for the removal of spurious variants based on the Kolmogorov-Smirnov statistical test. We experimentally validated our filters by comparing Pool-seq to individual sequencing data showing that the filters remove most of the false variants while retaining majority of true variants. The proposed guideline is fairly generic in nature and could be easily applied in other Pool-seq experiments.

A Genomic, Transcriptomic and Proteomic Look at the GE2270 Producer Planobispora rosea, an Uncommon Actinomycete.

We report the genome sequence of Planobispora rosea ATCC 53733, a mycelium-forming soil-dweller belonging to one of the lesser studied genera of Actinobacteria and producing the thiopeptide GE2270. The P. rosea genome presents considerable convergence in gene organization and function with other members in the family Streptosporangiaceae, with a significant number (44%) of shared orthologs. Patterns of gene expression in P. rosea cultures during exponential and stationary phase have been analyzed using whole transcriptome shotgun sequencing and by proteome analysis. Among the differentially abundant proteins, those involved in protein metabolism are particularly represented, including the GE2270-insensitive EF-Tu. Two proteins from the pbt cluster, directing GE2270 biosynthesis, slightly increase their abundance values over time. While GE2270 production starts during the exponential phase, most pbt genes, as analyzed by qRT-PCR, are down-regulated. The exception is represented by pbtA, encoding the precursor peptide of the ribosomally synthesized GE2270, whose expression reached the highest level at the entry into stationary phase.

High-throughput genetic characterization of a cohort of Brugada syndrome patients.

Brugada syndrome (BrS) is an inherited cardiac arrhythmic disorder that can lead to sudden death, with a prevalence of 1:5000 in Caucasian population and affecting mainly male patients in their third to fourth decade of life. BrS is inherited as an autosomal dominant trait; however, to date genetic bases have been only partially understood. Indeed most mutations are located in the SCN5A gene, encoding the alpha-subunit of the Na(+) cardiac channel, but >70% BrS patients still remain genetically undiagnosed. Although 21 other genes have been associated with BrS susceptibility, their pathogenic role is still unclear. A recent next-generation sequencing study investigated the contribution of 45 arrhythmia susceptibility genes in BrS pathogenesis, observing a significant enrichment only for SCN5A. In our study, we evaluated the distribution of putative functional variants in a wider panel of 158 genes previously associated with arrhythmic and cardiac defects in a cohort of 91 SCN5A-negative BrS patients. In addition, to identify genes significantly enriched in BrS, we performed a mutation burden test by using as control dataset European individuals selected from the 1000Genomes project. We confirmed BrS genetic heterogeneity and identified new potential BrS candidates such as DSG2 and MYH7, suggesting a possible genetic overlap between different cardiac disorders.

A glycosylated, labionin-containing lanthipeptide with marked antinociceptive activity.

Among the growing family of ribosomally synthesized, post-translationally modified peptides, particularly intriguing are class III lanthipeptides containing the triamino acid labionin. In the course of a screening program aimed at finding bacterial cell wall inhibitors, we discovered a new lanthipeptide produced by an Actinoplanes sp. The molecule, designated NAI-112, consists of 22 amino acids and contains an N-terminal labionin and a C-terminal methyl-labionin. Unique among lanthipeptides, it carries a 6-deoxyhexose moiety N-linked to a tryptophan residue. Consistently, the corresponding gene cluster encodes, in addition to the LanKC enzyme characteristic of this lanthipeptide class, a glycosyl transferase. Despite possessing weak antibacterial activity, NAI-112 is effective in experimental models of nociceptive pain, reducing pain symptoms in mice in both the formalin and the chronic constriction injury tests. Thus, NAI-112 represents, after the labyrinthopeptins, the second example of a lanthipeptide effective against nociceptive pain.

Whole genome SNP genotyping and exome sequencing reveal novel genetic variants and putative causative genes in congenital hyperinsulinism.

Congenital hyperinsulinism of infancy (CHI) is a rare disorder characterized by severe hypoglycemia due to inappropriate insulin secretion. The genetic causes of CHI have been found in genes regulating insulin secretion from pancreatic ß-cells; recessive inactivating mutations in the ABCC8 and KCNJ11 genes represent the most common events. Despite the advances in understanding the molecular pathogenesis of CHI, specific genetic determinants in about 50 % of the CHI patients remain unknown, suggesting additional locus heterogeneity. In order to search for novel loci contributing to the pathogenesis of CHI, we combined a family-based association study, using the transmission disequilibrium test on 17 CHI patients lacking mutations in ABCC8/KCNJ11, with a whole-exome sequencing analysis performed on 10 probands. This strategy allowed the identification of the potential causative mutations in genes implicated in the regulation of insulin secretion such as transmembrane proteins (CACNA1A, KCNH6, KCNJ10, NOTCH2, RYR3, SCN8A, TRPV3, TRPC5), cytosolic (ACACB, CAMK2D, CDKAL1, GNAS, NOS2, PDE4C, PIK3R3) and mitochondrial enzymes (PC, SLC24A6), and in four genes (CSMD1, SLC37A3, SULF1, TLL1) suggested by TDT family-based association study. Moreover, the exome-sequencing approach resulted to be an efficient diagnostic tool for CHI, allowing the identification of mutations in three causative CHI genes (ABCC8, GLUD1, and HNF1A) in four out of 10 patients. Overall, the present study should be considered as a starting point to design further investigations: our results might indeed contribute to meta-analysis studies, aimed at the identification/confirmation of novel causative or modifier genes.